arXiv: 1502.00813vl [physics.flu-dyn] 3 Feb 2015 


Mechanical and statistical study of the laminar hole formation in 

transitional plane Couette flow 

Joran Rolland* ** 

Submitted to Eur. Phys. J B: 28/05/2014, accepted: 03/02/2015. 


Abstract 

This article is concerned with the numerical study and modelling of two aspects the formation of laminar holes 
in transitional turbulence of plane Couette flow (PCF). On the one hand, we consider quenches: sudden decreases 
of the Reynolds number R which force the formation of holes. The Reynolds number is decreased from featureless 
turbulence to the range of existence of the oblique laminar-turbulent bands [R g ;i?t]- The successive stages of the 
quench are studied by means of visualisations and measurements of kinetic energy and turbulent fraction. The 
behaviour of the kinetic energy is explained using a kinetic energy budget: it shows that viscosity causes quasi 
modal decay until lift-up equals it and creates a new balance. Moreover, the budget confirms that the physical 
mechanisms at play are independent of the way the quench is performed. On the other hand we consider the natural 
formation of laminar holes in the bands, near R g . The Direct Numerical simulations (DNS) show that holes in 
the turbulent bands provide a mechanism for the fragmented bands regime and orientation fluctuations near R g . 
Moreover the analysis of the fluctuations of kinetic energy toward low values demonstrates that the disappearance 
of turbulence in the bands can be described within the framework of Large Deviations. A Large Deviation function 
is extracted from the Probability Density Function of the kinetic energy 

Keywords: Shear turbulence - Transition to turbulence - Nonequilibrium Statistical physics and Large Deviations, 

PACS: 47.27.nb - 47.27.Cn 05.70.Ln 


1 Introduction 

Using approaches from mechanics and statistical physics, this article examines the formation of laminar holes in 
featureless turbulence or in the oblique band regime of transitional plane Couette flow (Fig. [U (a)). Like the other 
shear flows, PCF, the flow between two parallel moving plates (Fig. |T| (b)), is always linearly stable. As a consequence, 
turbulence can only be triggered by finite amplitude disturbances. The transition is controlled by the Reynolds number 
of the flow. Above a Reynolds number R t — 415 the turbulence is said to be featureless and invades the whole domain 
PP- If the Reynolds number is decreased below R t and kept above R g ~ 325, the flow displays the intriguing oblique 
bands regime (Fig. [Q (a)) [2J |3J 0]. The coexistence of laminar and turbulent flow is spatially regular and statistically 
steady. Around i? g , the flow displays complex dynamics of local relaminarisations [5], fragmented bands Hi, spot 
splitting (7j and growth or retraction of spots [IHHli- 

Understanding the formation of laminar holes (or laminar troughs), that is region of very low kinetic energy in 
otherwise turbulent flow, is a key to explain the laminar turbulent coexistence. For that matter, two types of laminar 
troughs can be distinguished. On the one hand, one finds the troughs appearing in featureless turbulent flow, after a 
sudden decrease of the Reynolds number termed quench mm- These holes lead to the statistically steady oblique 
bands regime. On the other hand, laminar troughs can appear naturally inside the bands, and may lead to the collapse 
of turbulence m HUD m- The first type of laminar hole formation is easier to capture, since it is bound to happen in 
a finite time whenever the Reynolds number is decreased below Rt . Assuming that the mechanisms are very similar in 
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Figure 1: (a): Colour levels of the norm of the departure to the baseflow, result of a DNS in a periodic domain of size 
L x x L z = 110 x 72, at R = 350 (See Tab. [Dp for settings), (b): Sketch of PCF, introducing the domain size and 
velocity of the plates. 


both cases, the study of quenches can bring insight to the understanding of the collapse of turbulence and vice versa. 
These two points of view provide an occasion to test two different approaches of modelling. 

At the scale of the band, or that of the coherent structure, energy budgets m can be derived from Navier-Stokes 
equations. The comparison of the budget to the data extracted from DNS can give simple arguments to explain the 
behaviour of the flow. A mechanism can be proposed after the identification of the weight of dissipation (viscosity) 
and production in the variation of the kinetic energy. Production is caused by a central process in sustainment of 
turbulence: lift-up, that is, the extraction of energy from the mean flow by streamwise vortices [l7| (See |TB] for a 
recent review in Newtonnian and non-Newtonnian flows). 

An new point of view has been proposed to explain the super exponential dependence of the life time of turbulence 
in pipe flow: extreme fluctuations of the velocity fields leading to the collapse of turbulence M- The study of such 
rare events is often undertaken in the framework of Large Deviations (15] . This point of view completes the classical 
approaches taken to understand the disappearance of turbulence near R g : dynamical systems theory applied to the 
chaotic evolution of the self sustained process of turbulence mm or application of the spatiotemporal intermittency 
approach of turbulent convection mmmm- Indeed, it allows one to compute escape probabilities or lifetimes from 
the properties of the phase space of Navier-Stokes equations. In that spirit, long exponential tails of the Probability 
Density Function (PDF) of the kinetic energy have been measured in modelling of PCF [111 , and their asymmetry 
has been quantified [12]. Using an analogy with few degrees of freedom models, one can argue that this asymmetry 
is the precursor of a transition. However, several questions remain open. Does the collapse of turbulence in wall 
bounded flows fit in the general framework of Large Deviations ? Can all the concepts and methods associated to 
Large Deviations be used to model and predict the transition ? The study formation of holes near R g is adapted to 
answer these questions. 

The investigation of the properties of hole formation from these two points of views structures the organisation of 
the paper. The DNS procedures and the quantities used to investigate the flow are described first (§ H. The quench 
numerical experiments are then presented (§ 13), seen at the scale of the band (§ 13.21) and at the scale of the velocity 
streaks and the streamwise vortices (§GL3- We then suggest and validate a model describing the behaviour of the flow 
during the quench (§ H. We eventually move to the study of laminar holes appearing in bands near R g (§ [5]). These 
results are discussed together in section [G] 

2 Procedure 

2.1 Notations and numerical procedure 

Plane Couette flow is the flow between two parallel plates, at fixed positions y = ±h. The two plates move in opposite 
directions at velocities ±Ue x (Fig. |T| (b)). Making the system dimensionless simplifies comparisons with other flows, 
particularly with laboratory experiments: it largely reduces the number of control parameters. Moreover, provided the 
non-dimensionalisation is based on the proper physical criteria, the range of control parameter at which the transition 
occur is the same in all wall bounded flows [3]. Lengths are in units of h , velocities in units of U (the velocity of the 
top plate) and times in units of h/U. The dimensionless velocity field is ye x + v with v the departure from the laminar 
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Table 1: Table summarising the settings number (jj), the size used and the corresponding resolutions and thresholds. 

base flow. The Reynolds number R = hU/v, with v the kinematic viscosity, is the control parameter of the system. 
Together with the sizes L x and L z , it controls the behaviour of the turbulent flow. 

The spatial discretisation of the velocity field v and its integration in time is based on the code ChannelFlow 
written by John Gibson |23j . Fourier decomposition is used in the streamwise (x) and spanwise (z) directions, with 
N X)Z dealiased modes. Meanwhile, N y Chebychev polynomials are used in the wall-normal direction (y). In the 
plane, boundary conditions are periodic, while no slip conditions are used at the walls. An in plane resolution of 
N x ,z/L X ' Z = 4 and a wall-normal resolution of N y = 27 was found to be sufficient to perform a DNS of the flow at 
these Reynolds numbers [22 [23]. The code was used with a reduced resolution N x>z /L XtZ = 8/3, N y = 15 as a low 
order model of the flow. It is used when long simulations or repetitive statistics have to be performed (cumulated 
durations of several 100000/i/C/) [2 :23- The oblique bands regime was found to be well rendered at this resolution. 
The reduction of resolution affects mostly the transition thresholds. They are reduced from [i? g ;i? t ] — [325; 415], at 
N y = 27, to [i? g ; R t ] — [275; 355], at N y = 15 [22]• For each domain size used in this article, the resolutions used and 
the corresponding thresholds are summarised in table [Q 

2.2 Measurements 

In order to monitor the state of the flow, we compute the spatially averaged kinetic energy, e = \ J dxdydz || v || 2 , 
from the square of the departure from the laminar plane Couette flow. Besides, we compute the turbulent fraction 
/ (sometimes called the intermittency factor), using a method introduced in a earlier studies of the bands [21 Si .25] • 
This method is based on the identification the turbulent regions of the flow. The domain is divided in small boxes of 
size l x x l y x l z = 2 x 1 x 2, in which || v || 2 is locally averaged. The energy in each cell is compared to a threshold 
c = 0.0125 : if the local average of || v || 2 is above c, the cell is turbulent, otherwise, it is laminar. This threshold has 
been chosen in a former study [2j. The turbulent fraction is then defined as the fraction of cells which are turbulent, 
i.e. the relative volume occupied by turbulence in the flow. 

A quantity indicating whether the turbulence is organised in bands and which orientation do the bands take is 
necessary. For that matter, we take advantage of the sinusoidal modulation of turbulence of the band [2 El [6]. Once 
the fundamental mode of the bands (k x , ±fc z ) is identified, two similar strategies can be followed to define a so-called 
order parameter. The first one is adapted if the size of the system (L X:Z ) is large relatively to the wavelengths of the 
bands (2n/k XtZ ). One can then compute two Hilbert transforms of the signal measured around modes (k x , ±fc z ) [2 [5]. 
In numerical studies, the signal is the velocity field, while in the experimental studies, the signal is the light intensity. 
This yields two complex functions a + (x, z, t) and a_(x, z, t) which vary slowly in space relatively to the wavelength of 
the band. Their moduli give the amplitude of the modulation of the respective ± orientation at position (x, z) at time 
t while their phase give the relative shift of each patch of band. 

The second strategy is adapted if the size of the system is comparable to the wavelengths of the bands. Following 
Tuckerman & Barkley [3], we proposed to use the Fourier transform of the x component of the velocity field to define 
the order parameters a± = m±(t)e l ‘ t>± ^ by computing [2: 

1 f y=1 

m 2 ± = - / \u x \ 2 (k x ,y,±k z )dy , c/)± = &rg(u x (k x ,0,±k z ) . (1) 

1 Jy =-1 

The phases (f>± of the order parameters give the relative position of the respective orientation of the band in the domain. 
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Figure 2: Colour levels of the norm of the departure to the laminar baseflow || v || 2 during a quench experiment, in a 
domain of size L x x L z = 110 x 32, R = 370. (a): T = 0 , (b): T = 50, (c): T = 150, (d): T = 500. 


Meanwhile, the moduli of the order parameters give the amplitude of the modulation of the respective orientation. In 
this article, we follow this approach. Since we do not study the behaviour of the phases, we only use m±, termed the 
order parameters by an abuse of language. 

The order parameters are the ideal quantities to study the behaviour of the band and its appearance from featureless 
turbulence at R t . Unlike the turbulent fraction, they capture the spatial organisation of turbulence. This is particularly 
interesting when the two orientations of the band are competing. The presence of the system in one state or another 
can be easily determined by the ratio to+/to_: if it is large, the orientation + is observed, if it is small, the orientation 
— is observed. If it is of order 1, a mix of both orientations, a regular lozenge [9] or a more defective coexistence of 
laminar and turbulent flow [25], is observed. The order parameters may not be as efficient as the turbulent fraction 
to characterise the change of regime at R g [7]. However, they still yield a precise information as to the organisation of 
the flow in the so-called fragmented bands regime [4] . 

3 Quenches 

3.1 Principle 

The first and most simple way to obtain and study laminar holes appearing in turbulent flows is to perform so-called 
quench experiments [lQl ■ A featureless turbulent flow, at Rq > Rt ( Rq = 500 at N y = 27) is taken as an initial condition. 
It is obtained by integrating in time a random velocity field until it reaches statistically steady wall-turbulence. At the 
initial instant of the quench experiment, the Reynolds number is suddenly decreased to the so-called arrival Reynolds 
number R\ > R g . Note that R± can be above Rt- The simulation is performed at this constant value of the Reynolds 
number. 

Decorrelated featureless turbulent initial conditions are obtained by integrating in time a featureless turbulent 
velocity field over a duration of 300 h/U. This is useful for the study of the effect of the initial local organisation on 
the spatially averaged behaviour of the flow. 

Practically speaking, the quenches correspond to two possible situations. If the Reynolds number is decreased, 
and the velocity fields are unchanged, they correspond to an increase by Rq/Ri = (Uh/vo) / (Uh/v\) = v\/v o of the 
kinematic viscosity v. Meanwhile, if the velocity is rescaled by Rq/R\ = (Uoh/v)/(U\h/v) after the decrease of R, this 
corresponds to a decrease of the velocity of the plates by R\/Rq = Ui/Uq. The first approach is studied extensively 
in this article. The second is used in several cases in order to show that both approaches are equivalent. Although 
the decrease of viscosity can only be achieved in experiments where a cooling system controls the temperature and 
therefore the viscosity, it is the most systematic one numerically speaking. 

Both implementations depart slightly from a laboratory experiment. Indeed, in a laboratory experiment, the 
viscosity or the velocity profile have to adapt to the new imposed boundary conditions. Both adaptations occur in 
a time scale respectively proportional to the inverse of the temperature diffusion coefficient or the Reynolds number. 
Since ChannelFlow solves the equations for the departure to the laminar plane Couette flow, this raises other 
implementation questions such as using an actual discontinuity at t = 0, causing numerical instability, or taking into 
account the inertia of an experimental apparatus. By testing both the reduction of U and the augmentation of v, we 
will argue in section [4] that our procedure captures the main physical mechanisms. 

In that matter, Poiseuille flows are probably the most adapted to the simulation of realistic quenches. Indeed, the 
Reynolds number is decreased through a decrease of the pressure gradient or the flow rate, which adjust through the 
whole flow at the speed of sound. The corresponding duration is much smaller than the diffusion time. 
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Figure 3: Times series for quenches (i?o = 500, R\ = 370), in two domains of size L x x L z = 56 x 48 (full line) and 
L x x L z = 110 x 32 (dashed line), including a simulation with initial velocity rescaled by Rq/Ri, red), (a) : average 
kinetic energy e, (b) : turbulent fraction /. 


3.2 Visualisations and global measurements 

3.2.1 The quench 

The quench is first studied for a given arrival Reynolds number R\ = 370 in two domains of sizes L x x L z = 110 x 32 
(with and without rescaling v by Rq/Ri) and L x x L, = 56 x 48 (Tab.[I]ftl and (j2). The first domain accommodates 
one wavelength of the band, of fundamental wavenumbers k x = 27 t/ 110 , k z = 27 t/ 32 . The second one only displays an 
intermittent laminar-turbulent coexistence [25] ■ In order to test the effect of the initial organisation of the flow on the 
quench, five decorrelated uniformly turbulent initial conditions are used for both sizes. We consider the time series of 
kinetic energy (Fig. [3](a)) and turbulent fraction (Fig. [3](b)), and compare them to visualisations (Fig. [2]). 

We start from a turbulent initial condition (Fig. [2] (a)). The first stage of the quench takes place for times 0 < t < 80 
and corresponds to a sudden decrease of both e and /. During this stage, no laminar hole appears (Fig. [2] (b)), however 
the turbulent structures increase in size and decrease in intensity. The most striking feature of this stage is that data 
from the five numerical experiments collapse on a master curve (Fig. [3]). This indicates an apparent deterministic 
behaviour, independent of the initial organisation. This feature can be used to make a study as function of i?i of 
this stage: only one numerical experiment is necessary to capture the whole quantitative behaviour of the flow. The 
other interesting feature of this stage is that the behaviour of the flow does not depend on whether the system can 
accommodate bands or not. This suggests that this stage is dominated by a local process. Two types of behaviours 
appear: the average energy has an exponential-like decrease (Fig. [3] (a)), while the decrease of the turbulent fraction 
is linear (Fig. [3] (b)). This difference will be discussed in section [3731 

The second stage of the quench takes place for times 80 < t < 350. It corresponds to the formation of the laminar 
troughs (Fig. [2] (c)). The organisation of the troughs is apparently random and bears no trace of the turbulent band yet. 
This stage can be divided into two parts: first, both e and / undershoot. Then they grow back to the neighbourhood 
of their average value at this Reynolds number. Once they have reached this neighbourhood, they fluctuate around 
it and do not vary significantly. Unlike the first stage of the quench, the evolution is more random and there is no 
correlation between each experiment. The formation of laminar trough is discussed in more details in section [5] 

The last stage leaves no clear trace in e and / [25]. This stage corresponds to the reorganisation of the holes into 
the oblique laminar-turbulent bands (Fig. [2] (d)). In all three stages, the rescaling of the velocity field by Rq/R± only 
causes small decreases of the time scales and shifts of the minimum values of e and /. 

The duration of the last stage depends on the domain size L x x L z of the system and the complexity of the 
defects that can appear in the band. Its duration is approximately of 1000 h/U if the system accommodates one band. 
It is much longer if the system accommodates many bands: There is a competition between domains of + and — 
orientations |24j . This case shares many similarities with coarsening dynamics [271 . If the comparison is quantitative, 
one would expect scaling laws for the duration of the reorganisation stage as a function of the domain size L x x L z . 
The exponent depends on the type of system as well as its dimension. The comparison is all the more tempting that 
Ginzburg-Landau equations, the prototype equation of coarsening, appears to be a very good model of the measure 
of modulation a±(x , z, t) [4]. 
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Figure 4: Fitted slope a and growth rate r of the time series of (a) : turbulent fraction, (b) : average energy 


3.2.2 Effect of the arrival Reynolds number 

A study of the quenches for a large range of arrival Reynolds numbers, 340 < R± < 440 is now performed for both 
sizes, using one numerical experiment per size and Reynolds number. The same initial condition is used for each 
Reynolds number. Data sampled in quenches using a rescaled initial condition are added. Following the observation 
of the former section, the time series of the energy are fitted by a exp(—t/r) and those of the turbulent fraction by 
at + b. The parameters r and —1/a correspond to the characteristic decay times of kinetic energy and the turbulent 
fraction respectively. The parameters a and b corresponds to the values of e and / (respectively) in the featureless 
turbulence regime. One simply finds b < 1, as is expected from the featureless turbulent regime. The slope of the 
turbulent fraction a(R) increases toward zero with the Reynolds number (Fig. [I] (a)). It does not depend on the size of 
the domain in which the experiment is performed. The decay times t(R) of the kinetic energy increases very regularly 
with Ri (Fig. [4] (b)). A large incertitude appear at R\ = 365, for L x x L z = 56 x 48: a and r have to be averaged over 
5 experiments. This can be caused by unusually short first stages (Fig. [3] (a)). 

The change of size is found to shift the values of r by an additive constant. In both cases, for the lowest values of 
arrival Reynolds number Ri , the behaviour of a and r appears to be independent of the value of the Reynolds number 
of the initial condition R 0 . However, t and — 1/a diverge as R\ reaches Rq = 500. Indeed, for R\ close to i?o = 500, 
the flow changes very little when the quench is performed, leading to approximately constant values of e and /. Both 
—1/a and r show that rescaling the velocity field as the Reynolds number is decreased leads to shorter time scales, 
but the same behaviour as a function of R. 


3.3 Local measurements 


This section investigates the quench at small scale. It aims at determining fate of a small part of the flow, given its 
initial condition. For that matter, the domain is divided in small boxes of size l x x l z = 6 x 4 and thickness l y = 1 
(y > 0 or y < 0). The size of the cell is chosen so as to contain one or two velocity streaks. Around these values, 
no dependence of the result on the size was found. The state of the flow is evaluated in these boxes by averaging 
the kinetic energy in the box: this gives a set of averaged kinetic energies {e,} = ^ f hox || v || 2 drrdyd,?. Note that 
yb ej = e. The study is performed in a domain of size of L x x L z = 110 x 32 with an arrival Reynolds number of 
Ri = 370. The set of {ei(t)} is studied as a function of |ei(0)}, for t ranging from 0 to 250. 

The result is presented as clouds of data and reveals two phenomena (Fig. (5j(a)) . The initial line becomes an 
expanding cloud, and the cloud moves toward low values of e. The set of {e^} is not separated into two populations: 
one of smaller ej( 0) that would relaminarise and one of larger ej( 0) that would stay turbulent. This confirms the 
observation of adjustment (Fig. [2]). Note that during that stage, the distribution of {e^} is strongly positively skewed, 
with a relative skewness of approximately +1. 

Several measures are used to quantify this observation. The slope of the cloud a at time t is defined by least squares: 


miny^(ei(f) 

i 


aei( 0)) 2 => a = 


J2 t ei(t)ei( 0 ) 

E^O ) 2 


( 2 ) 
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Figure 5 : Local diagnostic of the quench R 0 = 500 —>• R\ = 370, L x x L, = 110 x 32, N y = 27. Cloud of points : (a) 
energy || v || 2 averaged in boxes at t = 0, t = 10, t = 50 as a function of that at t = 0. (b) : Slope of the clouds, (c) : 
correlation and dispersion coefficients of the clouds as a function of time. 


The normalised dispersion around the line is measured by the coefficients: 


£(*(*) - ^( O)) 2 


dispersion = i&e((o)2 + e((<)2 . 
And the correlation coefficient, i.e. the sense of a linear fit of the data, by: 

XXei(t) -e(t))(ei(0)-e(0)) 


correlation = 


VEfeW - e 0)) 2 EMo) - e(0)) 2 


(3) 

(4) 


The slope as a function of t is displayed in figure [5] (b). One can see the same exponential decay as that of the spatially 
averaged energy. The characteristic decay time is approximately 100, in agreement with data found at this Reynolds 
number (Fig. [4]). The cloud of {e.;} disperses and decorrelates in a similar duration f ~ 50 (Fig. [5] (c)), as can be 
visualised in figure [5] (a). This indicates that after this time, the slope measures the behaviour of the average energy 
rather than a local decay of each cell. 

Note that the two different behaviours of / and e in the first stage of the quench may be caused by the same 
phenomenon seen through different observables. Indeed, the kinetic energy e is governed by exponential decay at small 
and large scale. Meanwhile, the decrease of the turbulent fraction / is more gradual. It happens by small increments 
because the kinetic energy in a cell has suddenly decayed below the threshold c. 


4 Modeling the kinetic energy during the quenches 


The measurement of the kinetic energy: 

e = \ J (<E0 2 + K) 2 + H) 2 ) dxdydz 


(5) 


during the quenches can be compared to a budget derived from the Navier-Stokes equations. Using periodic boundary 
conditions and tensorial summation conventions with indices i,j, i' = {1, 2, 3} = {x, y, z}, the budget reads [15] : 


d t e = 


- / V X Vy 


^-dxdydz-^ J \djVi){djVi) J^di' . 


production 


dissipation, <0 


( 6 ) 


This equation has one production term, energy extraction from the base flow, and a viscous dissipation term. In 
the case of dimensionless plane Couette flow, the equation is simpler: U(y) = y and dU/dy = 1. When the flow is 
statistically steady, one approximately has dte — 0: the dissipation balances the production. 
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Figure 6: Production (blue lines) and dissipation rates (red lines) of the energy budget evaluated in a quench from 
R 0 = 500 to i?i = 370 in a domain of size L x x L z = 110 x 32, for an augmentation of viscosity (full lines) and a 
reduction of velocity (dashed lines). 


Since the two terms of the right hand side have a structure different from the left hand side, this equation is not 
tractable as such. One can either estimate both terms from DNS (Fig. [6]) or propose a simplified analytical treatment. 


We can first examine the behaviour of both term via DNS (Fig. (6]), using two decorrelated initial conditions, with 
and without rescaling the velocity field by Rq/Ri. When the quench is performed, 1 /R, and therefore the dissipation, 
is suddenly increased. As a consequence, one has c^e < 0. Both the dissipation and the production decay exponentially, 
until they once again balance one another at time t ~ 150. This corresponds to the end of the fast decrease of e, and 
to the slower readjustment to the new Reynolds number. In this stage, <9*e is much smaller (Fig. [jJJ and leads mostly 
to fluctuations. One can see that the behaviour of the production and the dissipation is the same for both quench 
procedures: in both cases, the exponential decay is driven by the sudden increase of dissipation between t = 0~ and 
t = 0 + . 

This behaviour can be described analytically. We first rewrite the budget. Indeed, the velocity v has periodic 
boundary conditions in x, z and vanishing boundary conditions in y. This allows one to decompose v on a base of 
Fourier modes in x , 0 and sine modes sin(((7rn)/2)(y +1)) in y , introducing k y = irn/2. Unlike Chebyshev polynomials, 
this choice is not adapted to the discretisation of the Navier-Stokes equations. However, this provides a tractable way 
of writing the energy budget. Indeed, in the case of plane Couette flow, one finds: 


e = 


\ : [{9j v i)(djVi) nd*' = ^( k Z + k l+ k x ,k v ,k. 


k x ,k z ,k x i 


k x ,k x ,k z i 


I dxdydz — ^ ' ^x,k x ,k y ,k z Vy,—k x ,k y ^k z 


(7) 


k x ,ky,k z 


We make the hypotheses: 

• The component v x dominates in the energy as well as in the dissipation (see Fig. [1}. 

• Instead of a wavepacket, the spectrum of v x is considered to be peaked around one given mode in k x , k z and k yi 
of wavelength corresponding approximately to the coherence length of the streak ~ 20, X z ~ 5. Note that 
the assumption of a single mode is not much affected by non-linear effects: indeed: the non-linearities conserve 
energy. The main drawback of this assumption is the possible dispersion of the wavepacket, caused mainly by 
the production term. 

The energy is then e = (l/2)\v x \ 2 (k x ,k y , kz), and the former equation reads: 




— \il. 7 ;~ — 


k l 


1 

2 


R 


( 8 ) 








Which can be rewritten and solved: 


kl + k 2 z + k 2 y ^ 


9t{\v x \) = -\v y \- 


|6x|(t) = |Dx|(0) exp — 


R 


+ j o \Vy\(t') ex P 


V ‘ R 

k 2 + k 2 + k 2 y 
R 


(t - 1’) At’ 


decay 


adjustment to v y 


(9) 


One finds exponential decay, caused by the initial unbalance between dissipation and production. This stops when 
extraction of energy from the base flow equals the viscosity. The adjustment to v y takes place over the same time scale. 
In practice, k 2 is negligible compared to k 2 (X z ~ 4) and k 2 . This yields a decay time proportional to I?: r = (1 /k 2 )R, 
with 1/k 2 of order one for both v x and e ~ v 2 . This is consistent with the measure of r for the smaller values of R 
(Fig. 0] (b)). The budget is invariant under a rescaling of velocity by Ri/Rq, corresponding to a quench by decrease 
of velocity. The quality of the approximation (shape of wavepackets, time evolution of v y ) is slightly impacted by this 
change. This explains the shift found in t(R) (Fig. [4] (b)). 

In the more realistic case of a time dependant profile U(y,t ) during the beginning of the quench, one has: 


\dt{ |w| 2 )-—|w| J v y ^(t)sm(k v (y + l))dy- 


* 2 , 


R 


( 10 ) 


The same treatment of the budget can be done, is replaced by / Ayv v ^(t) sm(k v {y + 1)) in equation ©. This 
leads to the same time evolution: exponential decay due to viscosity until lift-up equals it. This shows that apparently 
artificial quenches are a good model of a laboratory experiment, in the sense that they are driven by the same physical 
mechanisms. 

At a small scale, typically that of our averaging boxes, this description does not strictly hold. The boundary 
conditions are not periodic, and flux terms that cancel out in average will appear in the budget: 


dte+ lJ E dj(vjViVi) dxdz + - J d x (UviVi) di' + ^ J (1 — 8jj')vidyVi djd2 + J e?»(wjP) JJdi' 

/ i = r r. y. ■»/ i n , ='r r. 


J=X,Z 




momentum flux 


friction at boundaries 


incompressibility 


= - J v ^ v v~^ dxdydz - J {djv^djVi) JJdi'. (11) 


They correspond to advection of perturbation by the large scale flow, which usually excites turbulence jU [S] . The 
integrals corresponds to the sum over the cells introduced in section RHil They cause decorrelated and more random 
time evolution at small scale (Fig. [5]). However, the balance between viscous decay and extraction is expected to play 
the same role. One can view equation QT] as equation [6] with an additional zero average noise term. This noise term 
explains why the distribution of {e^} remains around a definite trend even though the dispersion is relatively large 
(Fig. [5] (a)). 


5 Holes in the bands 

We move from the formation of holes forced by the sudden change of Reynolds number to the natural apparition of 
holes inside the bands at a given Reynolds number. We present a case of change of orientation caused by the holes and 
reformation of bands. Then we propose a framework to study the fluctuations of quantities like the kinetic energy. 

5.1 Description 

An example displaying hole formation and of change of orientation at low Reynolds number is presented in figure [7] 
(colour levels of || v || 2 ) and figure [H] (time series). The domain studied has a size of L x x L z = 110 x 72 (See Tab. [TJ 
j)3). The band is prepared at R = 330, using a quench from R = 500 and letting it adjust for a duration of 2000. The 
Reynolds number is slowly decreased toward R = 317 in order to favour the formation of holes. 
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Figure 7: Norm of the departure to the laminar baseflow v 2 in the y = —0.62 plane at four successive times in a domain 
of size L x x L z = 110 x 72, at R = 317.. (a) : T = 3220, (b) : T = 3320, (c) : T = 3450. (d) : T = 3800 (see timescale 
on figure [8]). 


The domain contains one well formed wavelength of the band (Fig. [7] (a)). This corresponds to an order parameter 
dominating the other (Fig. [5] (a)) and a larger turbulent fraction (Fig. [8] (b)). After a long time (~ 3000 h/U), two 
holes appear at two different locations of the band (Fig. [7] (b)). This event is concomitant with a clear decrease of the 
turbulent fraction (Fig. [8] (b)) and a slight decrease of the dominating order parameter (Fig. [5] (a)). A small germ of 
spot is left isolated in the flow, which grows instead of decaying. The growth takes an oblique, asymmetric band-like 
form (Fig. [7] (c)), with an orientation opposite to that of the larger part of the band. This causes an increase of 
the turbulent fraction (Fig. [8] (b)), and a decrease of the dominating order parameter while the dominated one grows 
(Fig. [5] (a)). Eventually, this spot hits the rest of the band, which reorganises in the original orientation (Fig. [7] (d)), 
through the oblique advection of vortical structures by the large scale flow mm- 

These events can lead to changes of orientation near R g . They were seen in our low order modelling [25] as well 
as in experiments [3]. The same procedure of residency time measurements, based on the detection of the position 
of flow in phase space relatively to the two metastable states, can be used. The relevant observable is the pair of 
order parameters m±. However, the physical modelling (Eyring-Kramers Theory of mean first passage times [28] ) was 
specific and based on the fact that the system was well described by a gradient stochastic differential equation. The 
hole formation requires a more general approach. 

5.2 A general statistical approach ? 

The formation of laminar holes, and more generally, the excursion of the turbulent flow toward the laminar state is 
a very rare event for most Reynolds numbers [7] HU [T3 M- The rarity of these events is simply quantified by their 
probability (or their probability per unit time, rate of probability). From a physical point of view, these events usually 
correspond to a specific route followed by the flow 1111- 

Probability of rare events, long exponential tails, escape trajectories etc. can be computed in the framework of 
Large Deviations, a general approach of equilibrium and out of equilibrium statistical physics. It is centred around a 
Large Deviations principle (LDP). The LDP is simply the convergence of the logarithm of the pdf P toward a rate 
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Figure 8 : (a) 
figure [7] (L x x 



: Time series of the order parameter and (b) of the turbulent fraction centered around the events of 
L z = 110 x 72, R = 317). 





Figure 9: Processing of the pdf of the kinetic energy sampled with the low order modelling procedure (N y = 15), all 
functions are centred around the average of the energy, (a) : Logarithm of the pdf of the energy, at R = 310, for four 
different sizes, (b) : Logarithm of the pdf of the energy rescaled by the size L x x L z , at R = 310 for four different 
sizes, (c) : Rate function J of the energy, for three different Reynolds number, for the four different sizes. 


function J as a parameter 7 of the system goes to zero: 

J(x) = lim — 7 ln(P(a;)). (12) 

7— >0 

This is particularly interesting when the pdf is not a Gaussian. Indeed, when 7 is the inverse of the size or a number 
of particle, it goes beyond results like the central limit theorem. This theorem only leads to Gaussian distributions 
and provides a quadratic approximation of J around its minimum. 

The kinetic energy is a good example to demonstrate that this framework can be applied to PCF. The formation 
of a hole in the bands, which is the rare event, is obviously associated with an excursion of the average energy e 
toward low values. In PCF, the small parameter can be the inverse of the size of the system, and a rate function 
J = lim L x L z -yoo — 1 n(P)/(L X L Z ). In order to test that proposition, PDFs of e in domains of increasing sizes, going 
from L x x L z = 110 x 64 to L x x L z = 660 x 48 (Tab. HI |4 — 7), are sampled with the low order procedure. 

One can see that such a principle is well verified in the case of the pdf of the kinetic energy (Fig. [9] (a,b), where 
they are centred around (e)). The fluctuations of the kinetic energy decrease as the size increases. There is also a 
good collapse of In (P)/(L X L Z ) on a master curve, even for the smallest size. This curve is identified with the rate 
function. For a wide range of Reynolds number, the rate function is a parabola, which widens as the Reynolds number 
is decreased (Fig. [9] (b,c)). This means that fluctuations toward low or high values are equally probable. The rate 
function loses symmetry only for Reynolds number close to i? g , due to the changes in the phase space. This describes 
the exponential tail of the PDFs for small e, from which the probability of the formation of a laminar hole can be 
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extracted. In practice the loss of symmetry of the rate function shows to what extent the relaminarisation becomes 
probable as R is decreased toward R g . 

6 Discussion 

This article presented a numerical study of the formation of laminar holes in low Reynolds number turbulence of plane 
Couette flow. The hole formation was modelled using approaches from mechanics and statistical physics. 

The first part (§0 §SJ) was the use of so-called quenches, i.e. sudden reductions of the Reynolds number from R 0 
to R\. The DNS showed that the quantitative behaviour of the flow during the quench was rather independent of the 
way the quench was performed. Examination of the kinetic energy budget not only gave an explanation for the decay 
rate of the kinetic energy, but also confirmed that all ways of performing the quenches were equivalent. 

The second part (§ 0 of this study was the natural formation of laminar holes into the turbulent bands near the 
Reynolds number of disappearance of the bands. The holes appear when a local fluctuation is large enough to cross 
the basin boundary between laminar and turbulent flow and can cause a local relaminarisation. In this case dissipation 
is too strong to be balanced by the energy extraction m- The holes can be closed by the advection of vortices by 
the large scale flow mm or lead to a full relaminarisation of the flow HU. However, our DNS showed an interesting 
alternative scenario: if two laminar holes occur in the band at the same time, an isolated germ of band can grow with 
another orientation. This is a mechanism leading to orientation fluctuation different from what is found near R t , in 
which the turbulent fraction was nearly constant |25j . This is very likely the mechanism responsible for the fragmented 
band regime observed in experiments [4j. 

Besides, we verified that the extreme fluctuations toward low values of quantities like the kinetic energy could be 
placed in the specific framework of Large Deviations when the size of the domain diverges. This was done by showing 
the asymptotic convergence of the pdf in logarithmic form. This may provide general tools to compute average escape 
times, i.e. the lifetime of turbulence [14j . probabilities of relaminarisation. It may also provide approaches to determine 
paths between turbulent and laminar flow that are different from the deterministic methods [29( l30j . Indeed, they 
allow one to compute the full trajectory, including the crossing of the separatrix. A first stage of this study could 
consist in studying the stochastic models of laminar-turbulent coexistence m- Indeed, most of the analytical large 
deviations results are worked out from stochastic or Markov chain models. 
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